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Pseudogap and Superconducting Fluctuation in High-T c Cuprates: 
Theory beyond 1-loop Approximation 
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The pseudogap phenomena induced by the SC fluctuation are investigated in details. We 
perform a calculation beyond the 1-loop approximation. The SC fluctuation is microscopically 
derived on the basis of the repulsive Hubbard model. The vertex corrections are collected in 
the infinite order with use of the quasi-static approximation. The single-particle excitations, 
NMR 1/TiT, spin susceptibility and superconducting transition temperature are discussed. The 
important role of the vertex correction is pointed out for the single particle spectral function. 
On the other hand, the validity of the 1-loop order theory is confirmed for other quantities. 
We shed light on the essential nature of SC fluctuation leading to the pseudogap from the 
comparison with spin and charge fluctuations. 

KEYWORDS: Superconducting fluctuation; spin fluctuation; vertex correction; high-T c cuprates; pseudo- 
gap 



1. Introduction 

Since the discovery in 1986, 1 many anomalous aspects 
of high-T c superconductors have been clarified by enor- 
mous studies. Among them, the pseudogap phenomena 2 
(Fig. 1) have indicated a break down of the Fermi liq- 
uid theory, which is a universal view on the low energy 
excitations in metals. 3 The unusual nature has indicated 
an appearance of a new concept in the condensed matter 
physics. 4 

From the experimental point of view, the pseudo- 
gap was found in the magnetic excitation by the nu- 
clear magnetic resonance (NMR). 5 At present anoma- 
lies have been observed in many experimental probes 2 
which include NMR, 5-11 neutron scattering, 12 electric 



transport 
d ' troscopy 



13-18 



22-24 



optical spectrum, 19 21 tunneling spec- 
and angle resolved photo emission spec- 



troscopy (ARPES). 25-28 

Among many theoretical proposals, 29-35 we have 
adopted the "pairing scenario" in which the pseudogap 
is induced by the superconducting (SC) fluctuation. This 
scenario has been indicated by several experimental re- 
sults cited above. In particular, ARPES, which presents 
a momentum resolved information, has provided a strong 
circumstantial evidence. The pseudogap in the single par- 
ticle excitation first opens around the zone boundary k ~ 
(•7T, 0) and gradually extends to the whole Fermi surface 
as approaching to the superconducting transition. The 
amplitude of the excitation gap does not change through 
T c . The latter has been observed also in the tunneling 



spectroscopy 



and in the c-axis optical conductiv- 



ity. 19-21 These results have indicated a close relation be- 
tween the superconductivity and pseudogap. 

Since there are many approaches included in the pair- 
ing scenario, we have classified them into three kinds 
in order to avoid any confusion. 36 One is the BCS-BEC 
cross-over 37-40 formulated in the Nozieres and Schmitt- 
Rink theory. However, this approach is not relevant 
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for high-T c cuprates because this is basically justified 
in the low density system. We have adopted the sec- 
ond approach which is the perturbation theory with re- 
spect to the SC fluctuation. 36 ' 42-45 Then, the resonance 
scattering described by the electron self-energy plays an 
essential role. The 1-loop order theory which is called 
T-matrix approximation has been widely used and en- 
abled us to perform a non-phenomenological treatment 
on the microscopic Hamiltonian. This approach is com- 
plementary with the third one, phase fluctuation the- 
ory. The amplitude of order parameter is fixed in 
the phase fluctuation theory, which will be justified in 
the non-perturbative region with respect to the higher 
order mode coupling. 

It is our basic standpoint that non-Fermi liquid behav- 
iors in cuprates should be derived from the Fermi liquid 
state. Because cuprates behave as a conventional metal 
in the over-doped region and there is no discontinuity be- 
tween over-doped and under-doped regions, the origin of 
the anomalous behaviors should be inherent in the Fermi 
liquid state. Such understanding will be valuable for a 
universal understanding of the strongly correlated elec- 
tron systems including transition metals, organic super- 
conductors and heavy fermion systems. Then, the Fermi 
liquid theory will be a reasonable starting point. The 
microscopic study is suitable for this purpose since the 
characteristics of each system should be clarified from 
the microscopic point of view. 

In the previous studies, 52, 53 a microscopic theory for 
the SC fluctuation is developed on the basis of the re- 
pulsive Hubbard model. Then, the T-matrix and self- 
consistent T-matrix approximations for the SC fluctua- 
tion have been used in combination with the FLEX ap- 
proximation. 54 It has been shown that single-particle, 
magnetic and transport properties are explained in a 
coherent way. The doping dependence including the 
particle-hole asymmetry is also reproduced without using 
any phenomenological assumption. The electron-doped 
cuprates are basically weak-coupling superconductors 
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Fig. 1. The phase diagram of high-T c superconductors. The hor- 
izontal and vertical axes indicate the doping concentration and 
the temperature, respectively. "AF" , "SC" , and "PG" denote 
anti-ferromagnetic, superconducting and pseudogap state, re- 
spectively. The onset curve for the spin fluctuation (T=Tn) and 
that for the SC fluctuation (T=T*) show cross-over tempera- 
tures. 



and do not show pseudogap. The pseudogap state in the 
hole-doped region is described as a regime with spin and 
SC fluctuations, which have both cooperative and com- 
petitive aspects. We have shown that the pseudogap is 
essentially induced by the SC fluctuation. 

The purpose of this study is an examination of the 
theory based on the 1-loop approximation. The previous 
studies have attributed the pseudogap to the serious cor- 
rection arising from the 1-loop order term. 36 ' 42-45 ' 52 ' 53 
This fact involves a naive question from the theoretical 
point of view. Do the higher order terms play a more se- 
rious role ? There is no clear evidence for the validity of 
1-loop order theory because so-called Migdal theorem is 
not generally applicable. In this paper we present a the- 
ory beyond the 1-loop approximation in order to answer 
this question. This study will clarify the roles of higher 
order corrections. 

Theoretical studies beyond the 1-loop approximation 
have been developed mainly for the Lee-Rice-Anderson 
model 55 which is an one-dimensional model with charge 
fluctuation. Then, the electronic system coupled to the 
classical Gaussian field has been investigated. 56-60 The 
exact estimation is possible in the limit of infinite cor- 
relation length. 56 An approximation for the finite corre- 
lation length has been suggested by Sadovskii 56 and its 
validity has been investigated extensively. 57-60 The the- 
oretical method developed there has been applied to the 
two dimensional systems with charge fluctuation, 61 spin 
fluctuation 62 and SC fluctuation, 63 respectively. We will 
explain in §3 that this method is formulated in the mi- 
croscopic theory for the Hubbard-typc Hamiltonian. It 
will be discussed that the estimation based on the in- 
finite correlation length is applicable more robustly in 
the present case. We apply this method as well as the 
Sadovskii's method to the repulsive Hubbard model in 
§4. The attractive Hubbard model is also briefly adopted 
as another typical example in §5.1. 

We will show that higher order corrections are renor- 
malized to be small, even if higher order terms develop 
in the scheme of naive perturbation calculation. It is 



concluded that the 1-loop order theory is qualitatively 
valid in many aspects (§4.2-4). On the other hand, the 
importance of the vertex correction will be illuminated 
for the single-particle spectral function, for which the 1- 
loop and self-consistent 1-loop approximations provide 
qualitatively incompatible results. It is shown that the 
non-self-consistent 1-loop approximation is qualitatively 
correct for this quantity (§4.1 and §5.1). This is an ex- 
ample in which the partial summation included in the 
self-consistent theory provides an unphysical result. We 
will show that this conclusion is non-trivial by compar- 
ing with the case of spin fluctuation and that of charge 
fluctuation (§5.2). Indeed, the roles of vertex corrections 
are determined by the symmetry of the order parameter. 
We furthermore illustrate that the differences between 
the order parameters are qualitatively understood from 
the naive perturbation theory (§5.3). 

This paper is constructed in the following way. The 
1-loop order theory is shortly summarized in §2. Higher 
order terms are classified in §3. We introduce the terms 
focused in this paper and formulate the infinite-order 
theory. A justification for the adopted approximation is 
given in Appendix A and B. The roles of vertex cor- 
rections are investigated in §4. The single particle ex- 
citations (§4.1), superconducting T c (§4.2-3) and mag- 
netic properties (§4.4) are discussed in details. The ver- 
tex corrections beyond the formulation in §3 are briefly 
discussed in the estimation of T c (§4.3). The Sadovskii's 
method is used in §4.4, where the naive perturbation is 
performed within the 4-loop order and compared with 
the Sadovskii's method for a justification. The case of 
charge fluctuation and that of spin fluctuation are over- 
viewed in §5, where the essential nature of SC fluctuation 
is clarified. The summary and discussion are given in the 
last section §6. 

2. Theory in the 1-loop order 

In this section, the 1-loop order theory is summarized 
before discussing the vertex corrections. After the basic 
mechanism of the pseudogap was described on the basis 
of the attractive model, the microscopic theory has 
been developed on the basis of the repulsive Hubbard 
model. 52,53 ' 64 An approach for the infinite-?/ d-p model 
has been also developed. 65 Here we explain the formula- 
tion and typical results of the repulsive Hubbard model 
for a discussion in the following sections. 66 

The Hubbard Hamiltonian is expressed as, 

H = J2s(k)c^ a C ka + Uj2n iA n id . (1) 

We consider the square lattice and choose the following 
tight-binding dispersion e(k), 

e(k) = — 2i(cos k x + cos k y ) + At' cos fc x cos k y — fi. (2) 

The parameters are chosen as 2t — 1, t'/t — 0.25, U/t = 
4.2, and the doping concentration is chosen to be 5 = 
1 — n = 0.1, unless we mention explicitly. 

The SC fluctuation is generally described by the T- 
matrix which is the propagator of SC fluctuation. The T- 
matrix is represented by the ladder diagram as Fig. 2(a) 
where the effective interaction in the Cooper channel is 
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(a) 



(b) 
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Fig. 2. (a) The diagrammatic representation of the T-matrix. We 
denote the effective interaction between quasi-particles as V a - (b) 
The self-energy in the T-matrix approximation. 



represented by V a . In the attractive models, the first or- 
der perturbation has been used for the estimation of V a . 
42-45 j n orc j er t describe the superconductivity in the re- 
pulsive Hubbard model, we have to derive the effective in- 
teraction from the many body effects. Here we adopt the 
FLEX approximation 54 which has been widely used as a 
microscopic description of the nearly anti-ferromagnetic 
Fermi liquid. 67,68 The following calculation can be re- 
garded as an improvement of the FLEX approximation 
so as to take account of the SC fluctuation. Then, the 
effective interaction is described by the spin and charge 
susceptibility in the following way, 



V a (k,k') = U+-U 2 Xs (k-k') 



1 



U 2 Xc {k-k>), (3) 



where spin and charge susceptibility is estimated as, 
Xo(q) , x _ Xo(q) 



Xs(q) 
xo(q) 



l 



Xc(<?) 



u X o(qV 

-J2G(k + q)G(k). 



(4) 



(5) 



Here, the renormalized Green function is expressed as 

1 



G(k) 



(6) 



;„ - e(k) - £(fc) ' 

The self-energy in the FLEX approximation is obtained 
as S(/c) = Sp(fc), where 



E F (*) = £> n («)G(fc - «), 



(7) 



f 



(8) 



Vn(q) = U 2 [- X s(q) + ^Xc(q)-Xo(q)]- 

With use of these functions, the superconducting T c is es- 
timated by the following eigenvalue equation (linearized 
Eliashbcrg equation), 



A o 0(fc) = -^V/ a (fc,fc')|G(fc')| 2 0(fc')- 



(9) 



The transition temperature is determined by the crite- 
rion that the maximum eigenvalue is unity (A e =l) at 
T = T c . The FLEX approximation is regarded to be the 
mean field theory in the sense that the SC fluctuation 
is not taken into account. It should be stressed that the 
filiashberg equation is derived from the strong coupling 



theory for the superconductivity which is needed for the 
quantitative estimation of T c . This effort is important for 
the present purpose because the value of T c plays an es- 
sential role for the SC fluctuation. 36 Note that the weak 
coupling theory significantly over-estimates the value of 
T 

The T-matrix has been estimated by using the ex- 
tended Eliashberg equation as, 

T(k 1 ,k 2 :q) = cP(k 1 )t(q)(!>*(k2), (10) 
A(<?) = -J2 ^*( fc )K(fc - p)G( P )G(q - p)cf>(p), (12) 

k,p 

9=Y, ^ feMfe). (13) 



fei ,k 2 



where the wave function is normalized as 

5>(*)i 2 = 1 - 



(14) 



Since this estimation of the T-matrix is correct around 
9 = 0, the role of SC fluctuation is appropriately repre- 
sented. 52 It should be noted that the parameter A(0) is 
equivalent to the maximum eigenvalue of the Eliashberg 
equation, namely A(0) = A c . Thus, the divergence of the 
T-matrix is equivalent to the criterion in the Eliashberg 
equation. 

The self-energy correction due to the SC fluctuation is 
obtained within the 1-loop order as (Fig. 2(b)), 



Z s (k)=J2 T (k,k:q)G(q-k). 



(15) 



The total self-energy is obtained by adding it to the con- 
tribution from the spin fluctuation, 



E(fe) = E F (fc) + S s (fc). 



(16) 



The lowest order theory has been denoted as FLEX+T- 
matrix approximation (FTA) while the self-consistent 
calculation has been denoted as self-consistent FLEX+T- 
matrix approximation (SCFT). The basic role of the SC 
fluctuation is described by the FTA which is the 1-loop 
approximation with respect to the SC fluctuation. The 
SCFT is the self-consistent 1-loop order theory and fur- 
thermore includes the feedback effect on the magnetic 
fluctuation. The SCFT is needed to explain the trans- 
port phenomena in which the feedback effect plays an 
essential role. 53 

Figure. 3 shows the spectral function and the self- 
energy obtained by FTA. We clearly see the anoma- 
lous behaviors. The imaginary part of the self-energy 
(Fig. 3(a)) shows the scattering peak of |ImI7 R (fc,w) 
around u> = 0. This feature is qualitatively different from 
the conventional Fermi liquid where the self-energy be- 
haves as |ImX' R (fe, un)\ cx J 2 . In this sense, the Fermi 
liquid theory breaks down in the pseudogap state owing 
to the SC fluctuation. The scattering peak around u> = 
leads to the pseudogap around k— (ir,0) as is shown in 
Fig. 3(b). This is the basic mechanism of the pseudogap 
phenomena induced by the SC fluctuation. The pscudo- 
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gap gradually occurs with decreasing the temperature 
and/or doping concentration. 

Concerning the momentum dependence, the Fermi liq- 
uid behaviors appear around k= (tt/2,tt/2) which is de- 
noted as "cold spot" in contrast to the "hot spot" around 
k= (ir, 0). The w-linear behavior instead of a>-square be- 
havior is induced by the spin fluctuation. This momen- 
tum dependence induces the "Fermi arc" which has been 
observed in ARPES. 25-28 Note that these anomalous be- 
haviors are suppressed in the electron-doped cuprates 
where the SC fluctuation is very weak. 36 ' 52 
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have obtained the sufficiently large value of b in these 
regions. 52 On the contrary, the coherence length in the 
under-doped region has been estimated as £n = 3 ~ 5 in 
the unit of the lattice constant. This is the reason why 
the SC fluctuation plays an important role. 

At the last of this section, we comment on the differ- 
ence between FTA and SCFT. In the SCFT the dressed 
Green function is assigned in eqs. 5-15 instead of the 
Green function obtained by the FLEX approximation. 
Then, the pseudogap similarly appears in most of quan- 
tities, such as DOS, 52 although effects of the SC fluctua- 
tion are reduced. However, the situation is quite different 
for the single particle spectral function, which has been 
widely interested in the studies of high-T c cuprates. The 
SCFT does not show the pseudogap in this quantity on 
the Fermi surface, while it does slightly apart from the 
Fermi surface. This discrepancy will be discussed in §4.1 
and §5.1. This issue is particularly important because the 
pseudogap in the spectral function is a novel character- 
istics of the strong coupling superconductor as well as 
because it is an experimental evidence 25-28 for the pair- 
ing scenario. 

3. Vertex correction 

In this section we classify the vertex corrections and 
formulate the infinite order theory. Most part of this sec- 
tion is provided from the general point of view. The va- 
lidity of approximations adopted here will be examined 
on the basis of the microscopic theory explained in §2 and 
developed in §4, which is relevant for high-T c cuprates. 

3. 1 Classification of vertex corrections 



) s. (2) v ( 3a > . < 3b > 

<^r\ ^ x&s 



Fig. 3. (a) The imaginary part of the self-energy and (b) the spec- 
tral function on the Fermi surface. A, B and C show the results 
at k = (0.07tt, n), k = (0.15ir, 0.8tt) and k = (0.43tt, 0.45tt), re- 
spectively. The parameters arc chosen as U/t = 3.2, <5 = 0.1 and 
T = 1.2T C . 



Because the SC fluctuation is characterized by the 
T-matrix around q = Q = 0, the following expansion 
(TDGL expansion) is instructive for a qualitative under- 
standing. 

9 



t{q,Q) = 



t + bq 2 - (ax + ia 2 )fJ' 



(17) 



The parameter b is generally related to the coherence 
length £o as b oc £ 2 , which is estimated as £o oc wf/T c mf 
in the clean limit. Here T C MF is the transition tempera- 
ture in the mean field theory. We understand that suffi- 
ciently large value of b is obtained in the weak-coupling 
superconductor. Then, the total weight of the SC fluctu- 
ation is generally small. This is the reason why the fluc- 
tuation is negligible in conventional superconductors. In 
this sense, the over-doped and electron-doped cuprates 
are classified into the weak-coupling superconductor. We 



^ ^ ^ -®r 



Fig. 4. The series of vertex corrections focused in this paper. The 
solid and wavy lines represent the propagator of electron and SC 
fluctuation, respectively. 



We adopt a perturbation expansion with respect to the 
coupling between electrons and SC fluctuation, which 
is called loop expansion. In more general perturbation 
scheme for the Hubbard Hamiltonian, 69 there appear 
other terms which are not included in the loop expan- 
sion. However, we simply ignore them and focus on the 
singular terms in the vicinity of the phase transition. The 
vertex corrections to the Green function are mainly dis- 
cussed in this subsection. In addition to them, the vertex 
corrections to the pairing correlation function will be in- 
vestigated in §4.3. 

First of all, we show the Feynmann diagrams cor- 
responding to the self-energy investigated in this pa- 
per (Fig. 4). Fig. 4(1) is equivalent to Fig. 2(b), and 
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Figs. 4(2), (3a) and (4a) are included in the self- 
consistent T-matrix approximation. The other terms are 
classified into the vertex correction. The diagrams in- 
cluding the electron loop are not included in Fig. 4. It 
should be noticed that the terms like Fig. 5(a) should 
be included in the 1-loop order term owing to the defini- 
tion. Considering this fact, we understand that all of the 
ignored terms include the higher order mode couplings 
between SC fluctuations, which correspond to the non- 
linear terms in the Ginzburg Landau expansion (4th, 6th, 

8th order terms). For example, Fig. 5(b) and Fig. 5(c) 

include the fourth order vertex (Fig. 5(d)) and sixth or- 
der vertex (Fig. 5(e)), respectively. Therefore, the estima- 
tion of the terms in Fig. 4 is sufficient in the renormalized 
Gaussian region where the higher order mode coupling 
does not play any essential role, namely when the sys- 
tem is not too close to the critical point. We note that 
the mode coupling effect in the level of the self-consistent 
renormalization theory (SCR) 67 is included by itself if we 
determine the fluctuation propagator self-consistently as 
in §4.2-4. 

The diagrammatic expression in Fig. 4 is exact if we 
consider the Lee- Rice- Anderson model. Most of the pre- 
vious studies beyond the 1-loop approximation 56-63 have 
been devoted to this phenomenological model. The im- 
portance of these terms has been pointed out for the 
model with sufficiently long range attractive interaction 
where the calculation is reduced to the zero-dimensional 
problem. 70 We will estimate these terms as an approxi- 
mation for the microscopic model. 



(a) 



(b) 



(e) 



-ASS. 




(f) 



Fig. 5. Examples of diagrams which are not included in Fig. 4. 
(a) A diagram included in the 1-loop order term owing to the 
definition, (b), (c) and (f) are diagrams including higher order 
mode couplings, (d) and (e) are the forth order and sixth order 
vertex, respectively. 



The higher order mode couplings are surely important 
when the system is close to the critical point. However, 
considerable part of them is represented by the renor- 
malization of the fluctuation propagator. For example, 
the terms in Figs. 5(b) and (c) are the cases. These cor- 
rections are represented by the TDGL parameters (sec 
eq. 17) and the shift of T c . Therefore, the qualitative roles 
of SC fluctuation are not affected, although these correc- 
tions may be important for the estimation of T c . In this 
sense we have mentioned that the corrections in Fig. 4 
are important in the "renormalized" Gaussian region. In- 
deed, the renormalization of the fluctuation propagator 



has been investigated from 1980's motivated by the the- 
oretical issues on the BCS-BEC cross-over. 71 ' 72 Then, 
the higher order terms arising from the forth order ver- 
tex (Fig. 5(d)) have been investigated and concluded to 
be not important quantitatively. We will investigate the 
different kinds of mode coupling in order to discuss the 
doping dependence of T c in §4.3. They include the lower 
order terms with respect to the fluctuation propagator, 
rather than the terms investigated in Refs.71 and 72. 

There remain corrections including higher order mode 
couplings but not discussed above. We show the lowest 
order term in Fig. 5(f) which is the 4-loop order term. 
This term is more singular rather than the same order 
terms in Fig. 4. However, we will estimate this term, and 
conclude that this correction is small in the wide region 
of the pseudogap state (Appendix B). 

3.2 Quasi-static approximation 

According to these discussions, we consider the vertex 
corrections represented by Fig. 4. We adopt the quasi- 
static approximation in order to formulate an infinite- 
order theory. In this approximation, the thermal fluctua- 
tion is taken into account but the quantum aspect of the 
fluctuation is ignored. This procedure has been adopted 
in the studies of critical phenomena. In the formulation, 
the zero Matsubara frequency fi„ = of the fluctuation 
propagator is taken into account. This procedure is gen- 
erally justified at high temperature and/or in the vicinity 
of the critical point. 

Indeed, the character of SC fluctuation emarges in the 
applicability of the quasi-static approximation. The con- 
dition for the quasi-static approximation is generally de- 
scribed as uj sc <C T where T is the temperature and 
u> sc is the characteristic frequency of the SC fluctua- 
tion. Interestingly, u sc is expressed as iv sc ~ Te, where 
e = (T — T c )/T c is the reduced temperature. Then, the 
quasi-static approximation is justified when e <C 1. This 
fact should be contrasted to the case of AF spin fluctua- 
tion. Then, the characteristic frequency of the fluctuation 
is expressed as w s ~ E-ps s , where e s denotes the reduced 
temperature for the magnetic instability. It should be 
noticed that the factor appearing in co sc is smaller than 
that in lo s by the order T/E-p. This fact implies the better 
accuracy of quasi-static approximation for the SC fluctu- 
ation rather than that for the spin fluctuation. This will 
be confirmed by the microscopic calculation based on the 
repulsive Hubbard model (see Appendix A). While the 
quasi-static approximation for spin fluctuation is basi- 
cally justified at high temperature, that for SC fluctua- 
tion is effective in much lower temperature region. 

Another interpretation of the quasi-static approxima- 
tion is obtained from the application to the 1-loop order 
theory. In §2, we have shown the anomalous behavior of 
the self-energy induced by the SC fluctuation. It is easily 
confirmed that this anomaly arises from the quasi-static 
part of the T-matrix t(q, 0„ = 0), which is singular in the 
vicinity of T c . The remaining contribution from the dy- 
namical part basically induces the Fermi liquid behavior. 
Thus, anomalous contribution leading to the pseudogap 
is appropriately included in the quasi-static approxima- 
tion. Indeed, the applicability of quasi-static approxima- 
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tion is a requirement of the pseudogap phenomena be- 
cause the pseudogap in the single particle spectral func- 
tion appears when the former is more significant than 
the latter. Therefore, the quasi-static approximation is 
sufficient for the purpose of this paper, while the closing 
of pseudogap may be described inadequately. 

Owing to the quasi-static approximation, only the mo- 
mentum summation remains as a practical task. Al- 
though an exact estimation is still difficult, we can avoid 
this difficulty by considering the "dirty" region where the 
quasi-particle mean free path is much smaller than the 
GL correlation length, namely I <C £gl- In this region, 
the g-dependence of the Green function can be ignored, 
where q is the momentum of fluctuation. It should be 
noticed that quasi-particle in the "hot spot" has a re- 
markably short mean free path owing to strong scatter- 
ing due to the AF spin fluctuation. 73-75 Therefore, the 
condition I <C £gl is relevant in this region. Note that 
the superconductivity and magnetic properties are domi- 
nated by this region of the Fermi surface, while transport 
phenomena arc not. 53 

In the "dirty" region, the Green function in the n-loop 
order is equivalently estimated as, 

G (n) (k) = (-l)"A 2 ™(fc)G (0) (fc)™ +1 G (0) (-fc)™, (18) 

where A 2 (fc) = A 2 \c/)(k)\ 2 and A 2 = TJ2 q |*(«,0)|, which 
is proportional to the thermal weight of the SC fluc- 
tuation. The equivalence of n-loop order diagrams will 
be briefly discussed in §4.4 where a justification is ob- 
tained from the numerical estimation. We have defined 
the Green function without including the SC fluctua- 
tion as G(°\h). In the present case, the renormaliza- 
tion from the spin fluctuation is taken into account as 
G<°)(fc) = G F (k) = (kj n - e(k) - EfW)- 1 . 

Counting the topological factor n\, the summation of 
the infinite series is obtained as an Eular's error func- 
tion, 63 which is numerically tractable with use of the 
integral representation, 

G ^- G ^l \ + tA 2 Q GF (kw dL (19) 

In this procedure the properties of the SC fluctuation 
are represented by the single parameter A. We use this 
approximation in §4.1-3 where the single particle proper- 
ties around T = T c and the value of T c are investigated. 
A more sophisticated estimation will be provided in §4.4 
where the spatial fluctuation is taken into account. Note 
that the latter is reduced to eq. 19 in the limit £gl —* oo- 

4. Microscopic Theory in the infinite-loop order 

Hereafter, we investigate the role of vertex corrections 
on the basis of the repulsive Hubbard model. The basic 
formulation has been provided in §2 and §3. 

Before going on to the main issue, we summarize the 
advantages of the microscopic theory. (1) the properties 
of the SC fluctuation is microscopically derived. This 
advantage has enabled us to discuss the doping depen- 
dence without using any phenomenological assumption. 
(2) The renormalization of the quasi-particles is taken 
into account. Owing to the electron correlation, the en- 
ergy scale is renormalized to be smaller, and therefore the 
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pseudogap occurs around T ~ O.Olt ~ 100K. Note that 
the pseudogap in the attractive Hubbard model appears 
in much higher temperature region (see §5.1). (3) The d- 
wave symmetry of the superconductivity is naturally de- 
rived. The characteristic momentum dependence of the 
superconducting gap is derived, which is different from 
the simple form A(fc) oc cosfc x — cosk y especially in the 
electron-doped cuprates. 36 (4) The strong enhancement 
of the AF spin fluctuation is taken into account. The 
magnetic properties, which will be discussed in §4.4 are 
dominated by the spin correlation around q <~ (n, n). (5) 
The momentum dependence of quasi-particles is taken 
into account. The momentum dependent life time plays 
an essential role for the understanding of transport phe- 
nomena in the pseudogap state. 53 The remarkably short 
mean free path at "hot spot" is induced by the spin fluc- 
tuation, which provides a consistency with the assump- 
tion £ G l > I- 

In the following, we focus on the higher order correc- 
tions with respect to the SC fluctuation. The multiple 
scattering arising from the spin fluctuation has been in- 
vestigated motivated by the indication of Schrieffer. 76 
The previous studies 77,78 have shown that this correc- 
tion enhances the effective vertex, however does not af- 
fect qualitative roles of the spin fluctuation. We ignore 
the higher order coupling between the SC and spin fluc- 
tuations beyond the renormalization of Green function. 
The interests are focused on the roles of SC fluctuation. 

4-1 Single-particle properties 

In order to perform a infinite-loop calculation, we ten- 
tatively consider the neighborhood of T c and assume 
£gl 3> I- Then, the Green function is expressed as eq. 19. 
Correspondingly, the self-energy in the 1-loop and self- 
consistent 1-loop approximation is expressed as, 

Ss{k) = A(fc) 2 G F (-fc), (20) 

S s (k) = A(k) 2 G(-k), (21) 

respectively. We have ignored the frequency dependence 
of the coupling vertex <p(k), which gives no important 
effect in the vicinity of T = T c . In this subsection, we ig- 
nore the feedback effects on the SC and spin fluctuations 
in order to fix fluctuation propagators in each approxi- 
mation. A fully self-consistent calculation including the 
feedback effect is performed in §4.2-4. We have confirmed 
that qualitatively same results are obtained in the self- 
consistent calculation. 

Here, we take into account the inter-layer coupling 
in order to avoid the singularity in the exactly two- 
dimensional system. The TDGL expansion for the 
Lawrence-Doniach model is adopted as, 

t(9 ' 0) ^ + &g2 + 2 L(l-cos gz ) ' ^ 

with introducing the phenomenological parameter r. 
Here the TDGL parameters to and b are microscopi- 
cally determined. The phenomenological parameter is 
expressed by the anisotropy of the SC fluctuation as 
r = Cz/^O' wnere £z is the coherence length along the 
c-axis in the unit of the inter-layer spacing. The exactly 
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Fig. 6. (a) The spectral function at the hot spot, fc = (0.057r,7r). 
(b) DOS. The parameters are chosen as T = T C MF and r = 0.001. 



two dimensional and isotropic three dimensional systems 
correspond to r — and r = 1, respectively. The small 
value of r <C 1 is a underlying assumption of eq. 22. We 
choose r = 0.01 ~ 0.001 as a typical value. 

We show the results of the single particle spec- 
tral function in Fig. 6(a). The expressions 20, 21 and 
19 are denoted as FLEX+T-matrix, FLEX+SCT and 
FLEX+QSG, respectively. It these estimations, the effect 
of SC fluctuation is represented by the single parameter 
A. We find that the temperature dependence of A is not 
significant around T = T c owing to the inter-layer cou- 
pling. Therefore, the single particle properties are almost 
temperature independent, while they are rapidly altered 
below T c . 79 The experimental results have supported this 
result. 22 ^ 28 

It is shown that the gap structure in the spectral func- 
tion is suppressed in the FLEX+SCT, while it clearly 
appears in the non-self-consistent calculation (FLEX+T- 
matrix). This is the discrepancy explained in §2. We see 
that the non-self-consistent 1-loop approximation is qual- 
itatively correct. The calculation including the infinite 
order vertex corrections (FLEX+QSG) shows the pseu- 
dogap in the spectral function, although the gap struc- 
ture is reduced from the lowest order approximation. It 
is concluded that the self-consistent 1-loop order theory 
under-estimates the effects of SC fluctuation, while the 
non-self-consistent one over-estimates them. This qual- 
itative conclusion will be clearly understood from the 
perturbative point of view (§5.3). It may be expected 
that the self-consistent 1-loop approximation is always 
incorrect as is discussed for the single particle spectrum. 
In the following, we show that this expectation is invalid 



and the SCFT is quite appropriate for the macroscopic 
quantities. 

In the present case, the pseudogap in the infinite- 
loop calculation is broadened by the competition with 
the Fermi liquid behavior arising from the spin fluctu- 
ation. We note that the Fermi liquid behavior is basi- 
cally obtained in the nearly AF Fermi liquid theory, 52 
although some deviations from the conventional Fermi 
liquid are derived. 36 ' 67 As a result, the difference be- 
tween FLEX+SCT and FLEX+QSG is significantly re- 
duced. This is an origin of the unimportance of vertex 
corrections for the macroscopic quantities. 

The DOS is a typical one, which is shown in Fig. 6(b). 
It is shown that the pseudogap is similarly obtained in 
three calculations. We will show that the 1-loop approx- 
imation provides a rather different result if we assume 
s-wave superconductivity and neglect the spin fluctua- 
tion (§5.1). In the present case, the difference between 
approximations is significantly reduced by the <i-wave 
symmetry, which allows low-energy excitations, and by 
the spin fluctuation. Thus, better convergence of the loop 
expansion is expected in this microscopic theory rather 
than the phenomenological models. 

4-2 Phase diagram 

In this subsection, we estimate the doping dependence 
of T c . For this purpose we perform three kinds of calcu- 
lation. In the previous paper 52 we have performed the 
SCFT for the estimation of T c where eqs. 3-16 are self- 
consistently solved. Here we denote the calculation with 
cq. 21 instead of eq. 15 as the same abbreviation. In 
addition to that, we perform the self-consistent calcula- 
tion with eq. 19 instead of eq. 15 which is denoted as 
self-consistent FLEX+Gaussian approximation (SCFG). 
The calculation with eq. 20 instead of eq. 15 is de- 
noted as non-self-consistcnt FLEX+T-matrix approxi- 
mation (NSCFT). These calculations correspond to the 
self-consistent 1-loop, infinite-loop and 1-loop approxi- 
mation, respectively. In these calculations, the feedback 
effect on the spin fluctuation is self-consistently taken 
into account, which enhances the value of T c . 52 

In the previous paper, 52 the critical point has been de- 
termined by the criterion A c = 1 — e instead of A e = 1 as 
a phenomenological inclusion of the inter-layer coupling. 
In the present study, we explicitly take into account the 
inter-layer coupling with use of eq. 22. The value of r 
roughly corresponds to the phenomenological cut-off as 
e ~ Abr. Since the value of b is typically b ~ 5, the pre- 
viously used value e = 0.02 52 corresponds to r ~ 0.001, 
which is relevant for high-T c cuprates. Owing to the finite 
value of r, the logarithmic divergence in two-dimension 
is cut off and A at T — T c is obtained as 



A 2 = ^Arccos^ 



J^log-. 
47r6 r 



(23) 



Here q c is the cut off momentum along the a6-plane, 
which is chosen to be q c — ir. This procedure enables 
us to avoid the finite size effect which is serious around 
T = T C . 

First, we show the effect of inter-layer coupling in 
Fig. 7(a) where the results of SCFG are shown. Even 
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if the value of r differs by an order, the difference of T c 
does not exceed the factor 2 in the whole doping range. 
This is owing to the logarithmic dependence of A (see 
eq. 23). The qualitatively same doping dependence for 
r = 0.01 and r = 0.001 is an important property which 
allows a universal nature of the pseudogap 80 indepen- 
dent of the inter-layer structure. Although the value of 
r is remarkably small in BSCCO, 13 ' 81 this difference af- 
fects only quantitatively. 

It is shown that T c has a maximum around 5 ~ 0.11, 
while the critical temperature in the mean field theory 
T C MF develops with decreasing 6 in the under-doped re- 
gion (see Ref.36). This qualitatively different behavior 
between T c and T C MF is consistent with our understand- 
ing of the phase diagram. 36 That is, the onset tempera- 
ture of the pseudogap follows T C MF below which the SC 
fluctuation becomes active more and more. On the con- 
trary, the T c is suppressed in the under-doped region ow- 
ing to the fluctuation. We will show that the energy scale 
of the pseudogap is generally determined by T C MF and not 
by T c , independent of the inter-layer coupling (§5.1). 

It is predominantly believed that the inter-layer cou- 
pling is reduced with decreasing the doping. We have 
shown in Fig. 7(a) the result including the variation of 
r. Then, we have chosen the linear interpolation r = 
0.001 + 0.06(<5 - 0.08) so as to be r = 0.001 at <5 = 0.08 
and r = 0.01 at 5 = 0.23. We see that the value of op- 
timal doping is increased by the variation of inter-layer 
coupling. 

We see from Fig. 7(b) that the higher order corrections 
do not change the qualitative behavior of T c . Three ap- 
proximations give similar results. In particular, the close 
results are obtained between SCFT and SCFG. There- 
fore, it is concluded that vertex corrections beyond SCFT 
has no important role for the estimation of T c . This is 
not owing to the details of microscopic theory, since quite 
the same conclusion is obtained for the attractive Hub- 
bard model (§5.1). Note that the feedback effect through 
the spin fluctuation furthermore reduces the differences. 
Since the pseudogap in the single particle excitations 
is under-estimated in the SCFT, the T c is higher than 
SCFG. However, it will be shown that the T c in the SCFT 
is still under-estimated compared with the calculation in- 
cluding the vertex correction on the pairing correlation 
function. Details are explained in the next subsection. 

4-3 Vertex correction to the pairing correlation function 
Thus far, we have investigated the vertex corrections 
on the Green function, which are shown in Fig. 4. Ac- 
tually, more careful discussion is necessary in calcu- 
lating the many body correlation function. The Maki- 
Thompson (MT) term and Aslamasov-Larkin (AL) term 
on the electric conductivity are typical ones which have 
been investigated since 1970's. These corrections on the 
magnetic properties have been investigated in Rcfs.82- 
84 where the unimportance of MT-term in the (i-wave 
case has been pointed out. The microscopic treatment 
on the electric conductivity and Hall coefficient has been 
performed on the basis of the SCFT approximation. 53 
Then, it has been shown that the dominant correction 
is the indirect feedback effect because of the character- 
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Fig. 7. (b) Results of SCFG for various value of r. The diamonds 
are obtained by taking account of the doping dependence of inter- 
layer coupling, (a) The doping dependence of T c . The abbrevi- 
ations, SCFT, SCFG, NSCFT have been explained in the text. 
The times shows the results including the vertex correction on 
the pairing correlation function. The details and discussion will 
be given in §4.3. 



istic momentum dependence of high-T c cuprates. Then, 
the unimportance of the MT and AL terms has been 
concluded in contrast to the knowledges for the weak 
coupling superconductors. The correction to the Nernst 
effect has been also investigated on the basis of the SCFT 
approximation 85 where the MT term plays an important 
role. 






+ 
+ 




Fig. 8. The diagrammatic representation for the vertex correc- 
tions on the irreducible pairing susceptibility. 



Here, we investigate the vertex corrections to the pair- 
ing correlation function which is equivalent to the T- 
matrix except for a finite constant. The pairing corre- 
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lation function determines the SC fluctuation, and fur- 
thermore the value of T c . Since we find that the role of 
SC fluctuation is robust as is expected in §3, we focus 
on the correction to T c . We investigate corrections to the 
irreducible susceptibility in the Cooper channel, which 
are shown in Fig. 8. These corrections are regarded as a 
contraction of the mode coupling higher than the sixth 
order. However, it is easily confirmed that these correc- 
tions include the lowest order term in the sense of loop 
expansion, rather than the terms studied previously. 71, 72 
In order to estimate the corrections in Fig. 8 as well 
as those in Fig. 4, we collect the all order terms by using 
the same approximation as in §4.1 and §4.2. Then, the 
g-dependence of Green functions is ignored. This simpli- 
fication enables us to obtain a closed expression, which 
is the Eliashberg equation modified as, 



(24) 



Correspondingly, the function \(q) is obtained in the fol- 
lowing way, 



A(«) 



4?{k)V^k-p)Q{p,q)<l>{p)- 



(25) 



Here, the functions P(k) and Q(p, q) are defined as. 



P(k) 



Q(p, q) 



\G F (k)\ 



2 + tA(k) 2 
G F (p)G F (q 



dt, 



(26) 



p) 



A(p)2|G F (p)| 2 -A(p. 

A(p) 2 



g r\G F (q-pW 



/ { 

Jo 



\G F (p)\- 2 + tA( P y 



■}e-*dt, (27) 



\G F {q-p)\-z + tA(p-qY- 
respectively. 

The result of T c has been shown in Fig. 7(b). We un- 
derstand that the T c in the present calculation is larger 
than that obtained in the SCFT. Thus, the vertex correc- 
tions represented by Fig. 8 are larger than those included 
in the renormalization of Green function (Fig. 4). 

Note that this conclusion is expected from the pertur- 
bative point of view. We see that the lowest order vertex 
correction included in Fig. 4 is in the 3-loop order. On 
the other hand, the 2-loop order term exists in the cor- 
rections represented by Fig. 8. It is understood from the 
simple estimation that this 2-loop order term enhances 
T c . Therefore, the result in this subsection is consistent 
with the perturbation theory where the lowest order term 
is believed to determine the qualitative tendency. 

We have confirmed that the TDGL parameter b in- 
creases owing to the correction in Fig. 8. Therefore, the 
SC fluctuation is reduced by these corrections, while 
the increase of T c enhances the thermal fluctuation. We 
have confirmed that this effect is not serious, quantita- 
tively. Therefore, we ignore the corrections in Fig. 8 in 
the following discussion about the magnetic properties. 
We close this subsection by noting that the qualitative 
feature of the doping dependence is still not altered, as 
shown in Fig. 7(b). 



4-4 Magnetic properties 

At the last of this section, we show the temperature 
dependence of NMR 1/T\T and uniform susceptibility 
X s (0). They are typical quantities showing the pseudo- 
gap. 5-11 In order to discuss the magnetic properties, 
we perform another calculation proposed by Sadovskii. 
While we have only to consider the situation £ql ^ I for 
the estimation of T c , the spatial fluctuation may be an 
essential nature above T c . Sadovskii's method has pro- 
vided a simple and reasonable estimation for the spatial 
fluctuation. A justification will be obtained by the com- 
parison with the numerical estimation within the 4-loop 
order. A Monte Carlo simulation is a method of use, as 
has been performed for the phenomcnological model with 
spin fluctuation and that with charge fluctuation. 86 How- 
ever, we carry out an approximate but analytic treatment 
in order to obtain a clearer understanding and to com- 
pare several approximations in an equal footing. 

In general, a spatial fluctuation gives rise to the broad- 
ening of self-energy in each order. In the method pro- 
posed by Sadovskii, this effect is taken into account 
by linearizing the inverse of Green function as G F (k — 
q) = l/{\ijJ n — s(k) + v(k)q — T, F (k)) and replacing the 
qr-summation with adding an imaginary quantity. The 
imaginary quantity is chosen to be the inverse of typical 
length scale, namely Then, the self-energy in the 
n-loop order is expressed as, 



-(-iyx F ((-iyk)±ip(i)v(k)^lr 



(28) 



The integer p(i) is the number of SC fluctuation prop- 
agators enclosing the Green function appearing in the 
i-th order. The sign + (— ) should be chosen in the 
upper (lower) half plane. The quasi-particle velocity 
is obtained as v(k) = \v(k)\ where v(k) = v(k) + 
dReE F (k,mT)/dk. 

This method was first proposed for one-dimensional 
charge fluctuation systems as an exact solution. 56 Al- 
though this procedure has been proved to be not exact 
beyond 1-loop order, 57 a good accuracy as an approx- 
imation has been numerically confirmed in case of the 
complex order parameter. 58 ' 59 Here this method is ap- 
plied to the two-dimensional SC fluctuation. 

The infinite summation of eq. 28 is expressed as a re- 
cursive fraction, 



Es(*0 - 



(29) 



A(fe) 2 



-G F 1 (-k)+in(l)v(k)(: G l + - 



A(fc) 2 



G F 1 (fc)+m(2)«(fc)C G i + ... 

where n(i) = [(i + l)/2] and [...] is the Gauss symbol. It 
is easily confirmed that this expression is attributed to 
cq. 19 in the limit £gl — * oo. On the other hand, this 
expression is attributed to the zeroth order theory in the 
opposite limit £ql —> 0. Thus, the expression 29 can be 
regarded as an interpolation between the limiting cases. 

It should be noted that the applicability of this method 
to the two dimensional case will need more discussion. In 
the previous application for AF spin fluctuation 62 and 
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CDW fluctuation, 61 the displacement of the fluctuation 
propagator { -s +g2 c -ffg2 g-fig^ is imphcitly Per- 
formed. This displacement is accurate in the interested 
region q ~ 1/£gl- Owing to this procedure, the mo- 
mentum summation is attributed to the one-dimensional 
problem. Therefore, the accuracy of Sadovskii's method 
is similarly expected in the two-dimensional case. It 
should be stressed that the one-dimension is not particu- 
lar for the phenomenological Gaussian fluctuation model, 
while the quantum fluctuation is essential in the one- 
dimensional microscopic models. Note that two or higher 
dimensional models where the Sadovskii's method is ex- 
act have been proposed. 61 ' 70 
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Fig. 9. Results of the loop expansion, (a), (b), (c) and (d) show 
the imaginary part of the self-energy in the 1-, 2-, 3- and 4-loop 
order, respectively. The dashed lines denoted as S'sM show the 
results of Sadovskii's method. The indices (3a), (3b) and (4a-4d) 
correspond to the diagrams in Fig. 4. The parameters are chosen 
as U/t = 3.2 and T = 0.0084 where T = 1.035T C MF . 



In order to provide another ground, we perform the 
loop expansion up to the 4-loop order. Then, the quasi- 
static approximation is adopted and the dominator of 
Green function is linearized. The imaginary part of the 
self-energy in each order is estimated by the numeri- 
cal integration. The result of the Sadovskii's method 
(eq. 28) is shown together in Fig. 9. We see the quali- 
tatively same behaviors between the Sadovskii's method 
and naive loop-expansion in each order. 

Another interesting finding is that the crossing dia- 
gram and non-crossing one are almost equivalent in the 
results of loop expansion. For example, the difference be- 
tween Fig. 4(3a) and Fig. 4(3b) is almost invisible in 
Fig. 9(c). We see the same feature in the 4-loop order 
terms. This feature is obtained even if we do not lin- 
earize the dominator of Green function. 87 It should be 
stressed that this feature is an essential assumption of the 
Sadovskii's method. The results shown in Fig. 9 have con- 
firmed that this assumption is satisfied with remarkable 
accuracy. From these results we believe that Sadovskii's 



method is appropriate at least for qualitative discussions, 
although quantitative accuracy has not been proved. 

The numerical factor n(i) — [(i + l)/2] in eq. 29 is 
precisely obtained by counting the diagrams. This factor 
is replaced with n(i) = 1 in the self-consistent 1-loop ap- 
proximation (SCFT) and with n(i) = <5^i in the 1-loop 
approximation (NSCFT). The role of vertex corrections 
are clarified by the comparison between these calcula- 
tions. 
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Fig. 10. (a) The results of NMR 1/TiT obtained by several ap- 
proximations, (b) The uniform susceptibility. We denote the 
SCFT approximation fully including the spatial and dynamical 
fluctuations as SCFT (full). 



We show the results of NMR 1/TiT and uniform 
spin susceptibility in Fig. 10. Here the estimation of 
J2q<q l*(<7i 0)| i s numerically performed with the cut-off 
q c = ir/2. Because the inter-layer coupling is not impor- 
tant except for the vicinity T c , we perform the calculation 
in the two-dimension. The NMR 1/TiT is estimated by 
using the fluctuation-dissipation theorem, which is ex- 
pressed as 



= y2F ± (q)[-Im X f(q,L0 



(30) 



where 



F ±- (<?) = g { Al + 2B ( cos Qx + cos q y )} 2 + 
^-{A 2 + 2B(cosq x + cosq y )} 2 . 



(31) 
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The hyperfine coupling constants A\ , A2 and B are cho- 
sen as A 1 = 0.845 and A 2 = -45. 88 We have added 
the results obtained by the SCFT without using the 
quasi-static approximation. Then, the momentum and 
frequency summations are performed by using the fast 
Fourier transformation (FFT) and taking account of the 
same cut-off q c = n/2. 

It is clearly shown that the pseudogap phenomena ap- 
pear in all of the calculations in a similar way. NMR 
1/T]T shows a maximum above T c . On the other hand, 
the uniform susceptibility decreases with decreasing tem- 
perature from much higher temperature. This qualita- 
tively different behavior is one of the characteristics of 
under-doped high-T c cuprates. 2 ' 5,89 This is owing to the 
interplay of the AF spin fluctuation and SC fluctuation. 

We see that the vertex correction beyond SCFT is al- 
most negligible. The unimportance of the vertex correc- 
tion in the present case is very significant compared with 
the previous subsections. The result for \jT{F is partic- 
ularly surprising because this quantity is very sensitive 
to the approximation, owing to the strong enhancement 
of the spin fluctuation. This fact indicates an extremely 
precise cancellation of the vertex corrections for magnetic 
properties. 

It is noted that the same conclusion is obtained by 
neglecting the g-dependence of the Green function as is 
performed in §4.1-3. Then, the effect of spatial fluctua- 
tion is obviously neglected. This nature of approxima- 
tion should be contrasted with the Sadovskii's method 
where the spatial fluctuation is over-estimated. 63 It is 
expected that the unimportance of vertex corrections for 
the magnetic properties is a robust conclusion, since it is 
confirmed by two complementary approximations. 

At the last of this section, we comment on the closing 
of the pseudogap. As has been shown in §4.1, the pseu- 
dogap appears in the spectral function around T = T c . 
As the temperature increases, the gap structure is grad- 
ually broadened. In the present parameters, the pseudo- 
gap at the hot spot disappears around T = 0.005 which 
corresponds to the maximum of NMR \/T\T. We have 
confirmed that the closing of the pseudogap is yielded 
by the spatial fluctuation but not by the decrease of A. 
Indeed, the temperature dependence of A is weak since 
the increase of temperature enhances the thermal fluctu- 
ation. Therefore, the closing of the pseudogap occurs as 
a blurring but not as a shrinking. This is quite consistent 
with the experimental observation in ARPES. Note that 
the quasi-static approximation is not valid far above T c . 
Then, the dynamical fluctuation is also important as well 
as the spatial fluctuation. 

5. Phenomenological Theory: More General As- 
pects 

In this section, the comparison to spin fluctuation 
and charge fluctuation is performed on the basis of the 
phenomenological theory. We will see that the qualita- 
tive role of vertex corrections discussed in §4 is non- 
trivial and determined by the symmetry of order param- 
eter. An analysis of the infinite order calculation within 
the phenomenological theory has been reported for the 
charge fluctuation, 56 AF spin fluctuation 62 and SC fluc- 



tuation 



63,70 



respectively. Here, we provide a comprehen- 
sive view and shed light on the essential character of SC 
fluctuation. 



5. 1 In case of SC fluctuation 

First, we summarize the phenomenological theory for 
the SC fluctuation. The results on the single particle 
properties are clearly understood. In this section, we ig- 
nore the fc-dependence of A(fc) and phenomenologically 
assume the value of A for simplicity. This situation corre- 
sponds to the attractive Hubbard model which is briefly 
discussed at the last of this subsection. Since another 
scattering process is not taken into account here, the 
retarded Green function in the 1-loop, self-consistent 1- 
loop and infinite-loop calculations are described, respec- 
tively, 

w + e(k) 



Gf(*,w) 



(w + e(k))(u> - e(fc)) - A 2 
1 2 



e(fc)l 



1 



4A 2 
! -e(fc) 2 



U! - 



■e(k) 



(w + e(fc))(w-e(fc)) 
e(k) 



tA 2 



(32) 



(33) 



*dt{34) 



Strictly speaking, the condition I -C £ql is not relevant in 
the present case since I = 00. However, these expressions 
are sufficient for the qualitative discussion. 

The Green functions in eqs. 32 and 34 imply the rela- 
tion to the normal Green function in the superconduct- 
ing state. We see that eq. 32 is equivalent to the normal 
Green function in the BCS theory. The square of the or- 
der parameter in the latter is replaced by the thermal 
weight of the fluctuation. The expression in the infinite- 
loop calculation corresponds to the normal Green func- 
tion averaged by the Gaussian distribution of the order 
parameter. Quasi-particles propagate under the fluctu- 
ating order parameter which is approximated to be clas- 
sical and Gaussian in the present case. If the mean free 
path is shorter than the correlation length, the spatial 
fluctuation of the order parameter is not effective for 
quasi-particlcs. Therefore, a physically relevant result is 
obtained as eq. 34. Thus, the quasi-static approximation 
provides a clear interpretation. On the other hand, it is 
difficult to find a clear interpretation of self-consistent 
1-loop approximation (see eq. 33). This fact implies a 
potentiality that higher order terms are inappropriately 
taken into account in this approximation. 

This possibility has been realized in §4.1 where the 
single particle spectral function has been discussed. We 
illustrate the spectral function obtained from eqs. 32-34 
(Fig. 11(a)). The failure of the self-consistent 1-loop or- 
der theory is clearly understood. The single peak struc- 
ture appears in the self-consistent 1-loop approximation 
(SCT), while the double peak structure is clearly ob- 
served in the 1-loop (T-matrix) and infinite-loop (QSG) 
approximations. 63 The gap structure in the QSG is broad 
compared with the 1-loop approximation owing to the 
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Fig. 11. (a) The spectral function on the Fermi surface (e(k) = 0) 
obtained by cqs. 32-34. (b) The imaginary part of the self-energy. 
We have chosen the parameter as A = 0.2. 



Gaussian distribution of the pair field. Thus, the result 
obtained in §4.1 is qualitatively understood from this 
phcnomenological theory. 

This observation implies that the vertex corrections 
significantly cancel the higher order corrections included 
in the SCT. Note that the SCT applied to the attractive 
Hubbard model is one of the conserving approximations 
formulated by Baym and Kadanoff. The present calcu- 
lation is an example that the preservation of some con- 
servation low does not necessarily improve the accuracy 
of the calculation. The TPSC approximation has been 
proposed as an improvement where some conservation 
laws unsatisfied in the SCT are satisfied phenomenolog- 
ically. 90 ' 91 Then, the property of fluctuations and the 
coupling vertex to them arc basically discussed. However, 
the present diagrammatic method indicates that the im- 
provement of SCT should be performed by including the 
multiple scattering like Fig. 4. 

The characteristics of each approximation will be clar- 
ified by showing the self-energy in Fig. 11(b). It is shown 
that the anomalous scattering peak around u> — al- 
most disappears in the SCT. This is the reason of the 
broad single peak in the spectral function. The scatter- 
ing peak remains in the infinite order calculation, while 
it is reduced from the 1-loop order theory. 

Here, we show the DOS in Fig. 12. We see that the SCT 
is appropriate with regard to the DOS. This is because 
the spectral function slightly apart from the Fermi sur- 
face is not so different between the SCT and QSG. On the 
other hand, the result of the T-matrix approximation is 
rather different from that of QSG. This discrepancy has 
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not been shown in the microscopic theory (Fig. 6(b)). 
Then, the 1-loop order theory has been much improved 
owing to the d-wave symmetry and to the scattering from 
spin fluctuation. 



"T 1 T 




CO 

Fig. 12. Results of the DOS. The dispersion relation eq. 2 is used 
with t' It = and n = 1. 



At the last of this subsection, we provide results for the 
attractive Hubbard model in Fig. 13. Then, the Green 
function is expressed as eqs. 32-34 and the value of A 
and the propagator of fluctuation are self-consistently 
determined. Here the inter-layer coupling is taken into 
account with use of eq. 22. Fig. 13(a) shows that the 
vertex correction beyond the SCT is not serious for the 
estimation of T c . This conclusion is in common to the 
microscopic theory in §4.2. 

The relation between the excitation gap A at T = T c 
and the transition temperature in the mean field the- 
ory T C MF is shown in Fig. 13(b). We clearly see the uni- 
versal behavior independent of the indicator of three- 
dimensionality r = (£ z /£o) 2 - T C MF and A has a universal 
relation, although the value of T c depends on the inter- 
layer coupling. This is an important finding for the whole 
understanding of cuprate superconductors. The result in 
Fig. 13(b) shows that the energy scale of the pseudogap 
is not determined by T c , but by T C MF , which is robust 
for the variation of the inter-layer structure. 

It should be noted that this relation is derived in the 
renormalized Gaussian fluctuation region focused in this 
paper. Therefore, the experimental observation is not an 
evidence for the fixed amplitude expected in the phase 
fluctuation region. 46-51 

5. 2 Comparison to spin and charge fluctuations 

Next, the possibility of the pseudogap phenomena in- 
duced by spin fluctuation or charge fluctuation is dis- 
cussed in order to extract the essential properties of SC 
fluctuation. The possibility on the AF spin fluctuation 
has been investigated extensively in the studies of high- 
ly cuprates. 32-34,62 Although we will stress the invalid- 
ity of the quasi-static approximation for the spin fluc- 
tuation (Appendix A), this approximation is used for a 
comparison. The role of charge fluctuation has been in- 
terested in the context of the stripe order. 92 In our opin- 
ion, the charge fluctuation is not important for a unified 
understanding of high-T c superconductors. However, the 
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Fig. 13. (a) T c in the attractive Hubbard model. We show the 
results of SCT at r = 0.01 and QSG at r = 0.01 and r = 
0.001, respectively. The thin solid line show T c in the mean field 
theory, (b) The relation between T C MF and A at T = T c in each 
calculation and parameter. 
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These expressions are interpreted as the Green function 
in the ordered state averaged by the Gaussian distri- 
bution. 62 The essential difference in the formula is the 
power of x in the integrand, which is determined by the 
symmetry of the order parameter. The charge fluctua- 
tion has the ^-symmetry, while the SC fluctuation and 
spin fluctuation has the U(l)- and S'[/(2)-symmetry, re- 
spectively. They are translated into the Ising, XY and 
Heisenberg symmetry, respectively. 



charge fluctuation is an important example in the present 
purpose because the role of vertex corrections is quali- 
tatively different. The difference between the charge and 
spin fluctuations has been discussed by the Monte Carlo 
simulation. 86 We show that the SC fluctuation is classi- 
fied into the same class as that of the spin fluctuation. 

Here, the discussion is restricted to the commensurate 
case where the ordering vector is (tt,tt). The infinite se- 
ries for the Green function is obtained in the same man- 
ner as in §3. The results are expressed as, 



G c (k) = J2 



(In- 1)!! 
2™ 



A 2 "G(°)(fc)" +1 G (0) (fc + Q)" 



(charge fluctuation), (36) 
Gs(k) - Yl (2n 2 t 1)!! A s 2 "G (0) (fc)" +1 GW(fc + Q) n 

n 

(spin fluctuation), (37) 

respectively. The estimation of the numerical factor is 
not easy in case of the spin fluctuation since not only the 
longitudinal mode but also the transverse one exist. We 
obtain the factor ^ 2 " 2 t 1 - )!! from the recurrence formula. 
Here we have defined the thermal weight of the charge 
and spin fluctuation as A 2 = 9cT^2 q Xc(Q + 9,0) and 
A s = 9s T Y, q Xs(Q + 9,0) where g c and g s are the ef- 
fective coupling constants. In the RPA or FLEX for the 
Hubbard model, these coupling constants are equivalent 
to U. Summing up the infinite series, the Green function 




Fig. 14. (a) The spectral function on the Fermi surface obtained 
by the infinite-loop calculation for charge, SC and spin fluctua- 
tions, respectively, (b) The imaginary part of the self-energy, (c) 
The DOS. The parameters are chosen as A = A c = A s = 0.2, 
t/t = and n = 1. We have shown the DOS in the non- 
interacting case for a comparison. 



The spectral function obtained from eqs. 38 and 40 
are shown in Fig. 14(a). We see that the pseudogap is 
not observed in case of the charge fluctuation, while it 
is observed in case of the spin fluctuation. This is easily 
understand from the fact that the function exp(— a; 2 ) has 
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the maximum at x = 0, while the functions xexp(— x 2 ) 
and x 2 exp(— x 2 ) have the maximum at |a;| > 0. This 
is the reason why the SC fluctuation and spin fluctua- 
tion are classified into the same class. Fig. 14(b) shows 
that the scattering peak leading to the pseudogap in the 
spectral function almost disappears in case of the charge 
fluctuation. 

It should be stressed that qualitatively the same result 
is obtained for three kinds of fluctuations, when we adopt 
the 1-loop or self-consistent 1-loop approximation. The 
differences obtained above are essentially owing to the 
role of vertex corrections. It is understood from the above 
discussion that the qualitative role of vertex corrections 
essentially depends on the symmetry of order parameter. 
In this sense, the results on the SC fluctuation are non- 
trivial. The [/(l)-symmetry of the order parameter plays 
an essential role. In case of the charge fluctuation, the 
appropriate result is obtained from the self-consistent 1- 
loop approximation rather than the non-self-consistent 
one. 

It is shown in Fig. 14(c) that the charge fluctuation in- 
duces the gap structure of DOS, which has been focused 
in the previous studies. 55, 56,58-61 Thus, the pseudogap 
in the quantities including the momentum summation is 
robust for the properties of order parameter as well as 
those of approximations. The suppression of DOS around 
lo = is more remarkable when the symmetry of the or- 
der parameter is higher. 

Before closing this section, we comment on the realistic 
problem for high-T c cuprates. Then, the pseudogap in- 
duced by the AF spin fluctuation is not clearly observed 
in the FLEX approximation. We consider that this is not 
owing to the nature of the FLEX approximation, which 
is the self-consistent 1-loop order theory, but owing to 
the dynamical fluctuation which is important as will be 
shown in Appendix A. 

5. 3 Perspective from the perturbation theory 

The qualitatively different role of the vertex correc- 
tions discussed in §5.2 is clarified from the perturbative 
point of view. Fig. 15 shows the diagrammatic represen- 
tation of the loop expansion in case of the spin or charge 
fluctuation. We have shown the diagrams within the 3- 
loop order. Note that Figs. 15(1), (2a), (3a) and (3b) 
are included in the self-consistent 1-loop approximation, 
which is the FLEX approximation 54 in the microscopic 
theory. The other terms are classified into the vertex cor- 
rection. 

In order to make a discussion clear, we introduce a 
common expression for each term. In case of the charge 
and spin fluctuation, we denote 

£ (2a) (fc) = ]T ViQiWto) 

91,92 

xG{k + q 1 )G{k + q 1 +q 2 )G{k + q 1 ) 1 (42) 

£< 2b )(fc) =Y J V{q 1 )V{q 2 ) 

91,92 

xG{k + q 1 )G{k + q 1 +q 2 )G{k + q 2 ) 1 (43) 



0) 



(2a) 




a) (3b) ( 3c ) ^_ X - N _ > ^ 



(3e) 



(3!) 



-n^S r^Vv 

W (3h) (3 i) 




Fig. 15. The diagrammatic representation of the self-energy 
within the 3-loop order. The term (2b) is the lowest order vertex 
correction in case of the charge or spin fluctuation. The term (3j) 
is the lowest order vertex correction in case of the SC fluctuation. 



£ (3a) (fc)= E V(q 1 )V(q 2 )V(q 3 )G(k + q 1 ) 

91,92,93 

xG(k + qi + q 2 )G(k + qi + q 2 + q 3 ) 
xG(k + qi +q 2 )G(k + qi ), (44) 
E< 3 j)(fc)= £ V(q 1 )V(q 2 )V(q 3 )G(k + q 1 ) 



xG(k + qi + q 2 )G(k + qi + q 2 + q 3 ) 

xG(k + q 2 + q 3 )G(k + q 3 ), (45) 

and so on. The vertex V(q) is defined as V(q) = ^g 2 Xc(q) 
and V(q) = | g 2 Xs(q) for charge and spin fluctuation, re- 
spectively. We have chosen the notation in which the co- 
efficient of the terms included in the self-consistent 1-loop 
approximation is unity. It is easily confirmed that odd 
order terms generally introduce the anomalous contribu- 
tion leading to the pseudogap, while even order terms 
reduce it. 

We see that the lowest order vertex correction is in the 
second order (Fig. 15(2b)). In case of the charge fluctu- 
ation, all of the terms contribute to the self-energy with 
a coefficient 1. Then, the self-energy in the second order 
is £( 2 >(fc) = £< 2a )(fc) + Z (2b \k) and that in the third 
order is S< 3 )(fc) = £J Q=a i;( 3Q: ) (k). In this case, the lowest 
order vertex correction reduces the anomalous contribu- 
tion. Therefore, the effect of the charge fluctuation is 
still overestimated in the self-consistent 1-loop approxi- 
mation. 

On the other hand, the coefficients are complicated 
in case of the spin fluctuation. We obtain the second 
order term as £< 2 )(fc) = S< 2a )(fc) - \^ 2h \k). The neg- 
ative coefficient —1/3 arises from the coupling between 
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the longitudinal and transverse modes. Thus, the sign of 
the lowest order vertex correction is opposite to the case 
of charge fluctuation. Therefore, the qualitatively differ- 
ent role of the vertex correction is expected. The self- 
consistent 1-loop approximation underestimates the ef- 
fects of spin fluctuation. Note that the 1-loop approxima- 
tion overestimates them since the vertex correction term 
-i£( 2b ) (k) has smaller coefficient 1/3 than £< 2a ) (k). The 
Monte Carlo study has reported the qualitatively consis- 
tent result about the difference between the charge and 
spin fluctuations. 86 

The numerical factor —1/3 has played an important 
role in the study of the pairing interaction. 77, 78,93 The 
lowest order vertex correction enhances the pairing in- 
teraction mediated by the spin fluctuation, 77, 78 while it 
reduces the pairing interaction mediated by the charge 
fluctuation. 93 Thus, the spin fluctuation is likely to be ef- 
fective for the anisotropic superconductivity rather than 
the charge fluctuation, owing to the vertex corrections as 
well as the numerical factor 3. 

At last, we clarify the case of SC fluctuation. Then, the 
diagrams in Figs. 15(2b) and (3c-i) are absent. The num- 
ber of the vertex correction terms is remarkably small in 
this case. Indeed, this is a characteristic property of the 
complex order parameter. We introduce the expressions 
for the remaining terms as, 





charge 


SC 


spin 


2nd order 


1 





-1/3 


3nd order 


8 


1 


-22/9 



£ (2a) (*0= £ WW 



(46) 



xG(-k + qi)G(k -qi + q 2 )G{-k + gi), 

£ (3a) (*0 = E www 

xG(-k + qi )G(k - qi + q 2 )G(-k + qi-q 2 + q 3 ) 
xG(k-q 1 +q 2 )G(-k + q 1 ), (47) 

s (3j) w= £ WWW 

91,92,93 

xG(-k + qi )G(k -qi + q 2 )G{-k + qi - q 2 + q 3 ) 



xG(k + q 2 -q 3 )G(-k + q 3 ), 



(48) 



where V(q) = t(q,0). The second and third order terms 
are obtained as ^(k) = £( 2a )(fc) and £( 3 )(fc) = 
£( 3a )(fc) + £( 3b )(fc) + £( 3 J)(fc), respectively. Because the 
coefficient of each term is unity, the lowest order ver- 
tex correction 'E^(k) has qualitatively the same role 
as ^ 3a, \k) which enhances the pseudogap phenomena. 
Thus, the role of vertex corrections clarified in §4 is ex- 
pected from the perturbative point of view. Because the 
vertex correction in this case is higher order than that of 
spin and charge fluctuations, it is expected to be less im- 
portant. This is an underlying origin of the unimportance 
of the vertex correction on the macroscopic quantities. 

Summarizing, the results obtained in §5.1 and §5.2 are 
qualitatively consistent with the expectation from the 
perturbation theory. The qualitatively different role of 
vertex corrections is explained from the lowest order ver- 
tex correction, even if the contribution from the higher 
order term is larger. Thus, the perturbation theory pro- 
vides another understandings for the results, although 



Table I. The coefficient of vertex correction terms. The cases of 
the charge, SC and spin fluctuation are shown, respectively. We 
show the coefficient of the second order term E' 2b ) and the sum- 
mation of the coefficients of third order terms S' 3c ^-E( 3 ^. 



we have to perform an infinite order calculation in order 
to obtain the results in a closed form. 

6. Summary and Discussion 

In this paper, we have investigated the role of higher 
order corrections beyond the T-matrix approximation for 
SC fluctuation. The results obtained from the detailed 
analysis are summarized in the following way. The first 
one is the importance of the vertex correction for the sin- 
gle particle spectral function. Then, the self-consistent 
T-matrix approximation is qualitatively incorrect, while 
the non-self-consistent one is appropriate. The second 
one is the unimportance of the vertex correction for 
the macroscopic quantities. We have explicitly estimated 
the DOS, superconducting transition temperature, NMR 
\jT\T and magnetic susceptibility on the basis of the re- 
pulsive Hubbard model. It has been shown that the self- 
consistent 1-loop order theory is rather precise for these 
quantities. We have found that the cancellation of the 
vertex correction for the magnetic properties is especially 
precise. These results basically provide a justification of 
the 1-loop order theory, although a care is necessary for 
the interpretation of ARPES measurements. The present 
study provides a clear point of view to the previous stud- 
ies within the 1-loop order calculation. 36, 42 ~ 45, 52,53 

We have pointed out that the d-wave symmetry of the 
superconductivity and the renormahzation from the spin 
fluctuation significantly reduce the effect of vertex cor- 
rections. This fact is natural since the higher order terms 
are generally suppressed by these properties. Therefore, 
a better convergence of the loop expansion is expected 
in the realistic situation rather than in the situation 
adopted in the phenomenological models. 

It should be stressed again that the justification of the 
1-loop order theory is obtained even when the naive per- 
turbation calculation fails. Then, the higher order cor- 
rections are renormalized by summing up infinite series. 
Note that the cancellation occurs between the even order 
terms and odd order ones which have qualitatively oppo- 
site behaviors at low energy. The same order terms en- 
hance rather than cancel each other. Therefore, the cal- 
culation within the finite order is remarkably inaccurate. 
The results in 3-loop order, 4-loop order seriously os- 
cillate. The 1-loop order theory is significantly precise 
rather than the other finite order calculation. However, 
we have pointed out in §5.3 that the qualitative roles of 
vertex corrections are generally determined by the lowest 
order term. The unimportance of the vertex corrections 
arising from the SC fluctuation rather than that from 
the charge and spin fluctuations is also expected from 
the perturbative point of view. Thus, the perturbation 
theory is still valid in this qualitative sense. 
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Note that the 1-loop order theory has been used widely 
in the fluctuation theory including the FLEX approxima- 
tion. 36 ' 54,67 The simplicity of the 1-loop order theory has 
enabled us to investigate rich issues in a coherent way. At 
the moment, there is no conclusive evidence for the va- 
lidity because so-called Migdal theorem is not generally 
applicable. It may be considered that this problem is se- 
rious in the strongly correlated electron systems because 
the strong coupling nature of the fluctuation is a char- 
acteristic property of them. However, the present study 
implies a wide applicability of the 1-loop order theory, al- 
though it inevitably fails in the deeply critical region. For 
example, we believe that the 1-loop order theory for the 
spin fluctuation is more effective than expected from the 
lowest order estimation for the vertex correction. 62 ' 77 ' 78 
In many cases, the contribution which is not included in 
the fluctuation theory is more important rather than the 
vertex correction, as a quantitative correction. 36 We ex- 
pect that the observation in the present study will be a 
typical example in the rich field of the fluctuation theory 
in the correlated electron systems. 
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Appendix A: On the validity of quasi-static ap- 
proximation 

This Appendix is provided as a ground for the quasi- 
static approximation which is used in this paper. For a 
quantitative estimation needed for this purpose, we show 
the results obtained in the microscopic theory developed 
in §4.4. Then, the quasi-static approximation has been 
used in combination with Sadovskii's method. We show 
the result of SCFG at T = 0.0035. Fig. A-l shows the fin- 
dependence of the pairing correlation function at q = 0. 
The spin correlation function at q — (ir, w) is shown for 
a comparison. In order to perform a comparison in an 
equal footing, we show the functions \t{q)/g\ — \\{q)/{l — 
\(q))\ and Ux s {q) = Uxo(q)/(l - Uxo(q)) and choose 
the parameters so that the values at £!„ = arc almost 
equivalent. 

It is clearly shown that the pairing correlation func- 
tion rapidly decreases for a finite Matsubara frequency. 
Therefore, it is expected that the dynamical part of the 
T-matrix does not play any important role. This conclu- 
sion has been confirmed also in the explicit calculation 
within the 1-loop order theory (§4.4). Then, the SCFT 
approximations with and without quasi-static approxi- 
mation provide similar results for magnetic properties, 
we have confirmed this fact also for the single parti- 
cle spectral function, although the result has not been 
shown. 

On the other hand, the spin correlation function takes 
not so small value at finite Matsubara frequency. In other 
words, the characteristic frequency of the spin fluctua- 
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Fig. A-l. The Matsubara frequency dependence of the T-matrix 
at q = and that of the spin susceptibility at q = (tt,tt). The 
functions \t(q)/g\ and Uxs(q) are shown for an equitable com- 
parison. The SCFG is adopted with T = 0.0035. 



tion remains to be larger than the temperature. Thus, 
it is concluded that the quasi-static approximation for 
the spin fluctuation is not valid for the relevant region in 
the pseudogap state. We wish to stress again that such 
a remarkable difference between the SC fluctuation and 
spin fluctuation is not owing to the enhancement factors 
1 / ( 1 — A (q) ) and 1/(1 — Uxo (<?)), but owing to the essen- 
tial property of each fluctuation. As has been explained 
in §3.2, the characteristic frequency of the fluctuation is 
expressed as ui sc ~ Te and ui s ~ Eps s for SC fluctuation 
and for spin fluctuation, respectively. We see a smaller 
coefficient in case of the SC fluctuation. Indeed, this is 
one of the characteristics of the phase transition with 
logarithmic divergence at T — > 0. For example, the su- 
perconductivity and SDW with perfect nesting are the 
cases. 

It should be noted that the validity of the quasi-static 
approximation is related to the formation of the pseudo- 
gap in a straightforward manner. While the anomalous 
contribution leading to the pseudogap arises from the 
quasi-static part, the dynamical part basically smears it. 
The validity of the quasi-static approximation is actually 
an underlying origin of the clear pseudogap induced by 
the SC fluctuation. 

Note that the invalidity of the quasi-static approxima- 
tion for the spin fluctuation is concluded in the pseu- 
dogap state, but it is not a general consequence. This 
approximation for the spin fluctuation theory will be 
appropriate in higher temperature region and/or suffi- 
ciently close to the magnetic instability. In the present 
case, the pseudogap in the magnetic excitation itself is 
an origin of the invalidity, since it means the reduced 
damping of spin fluctuation. 

Appendix B: Self-energy represented by 
Fig. 5(f) ' 

In this paper, we have ignored the contribution from 
the non-Gaussian terms which include higher order mode 
couplings. Although considerable part of them is repre- 
sented by the renormalization of fluctuation propagator 
(see Figs. 5(b) and (c)), the role of remaining terms is 
not clear. The lowest order term included in the remain- 



J. Phys. Soc. Jpn. 



Full Paper 



Youichi Yanase 17 



10 
3 5 

lo 



1 1 1 1 

- 


1 1 1 


1 ---Fig.5(f) 

Fig.4(4a) " 




i U i 


i 



-0.02 



-0.01 





CO 



0.01 



0.02 



Fig. B-l. The imaginary part of the self-energy. The dashed line 
shows the contribution from the non-Gaussian term represented 
by Fig. 5(f). The solid line shows the self-energy represented by 
Fig. 4(4a). The temperature is chosen as T = 0.004. 



ing part has been shown in Fig. 5(f). Here, we provide a 
ground for our disregard of these terms. 

We explicitly estimate the self-energy Fig. 5(f) and 
compare with the 4-loop order term in Fig. 4. Here, 
we determine the T-matrix and Green function by us- 
ing the SCFG as is explained in §4.4. For the estimation 
of Fig. 5(f), we simply ignore the frequency and momen- 
tum dependences of the 4th order vertex T represented 
by Fig. 5(d), as is done in the Ginzburg-Landau theory. 
Since there are two equivalent diagrams of Fig. 5(f) ow- 
ing to the spin index, we multiplies the non-Gaussian 
terms by the factor 2. 

The result is shown in Fig. B-l. It is understood that 
the non-Gaussian term is much smaller than the same 
order vertex corrections taken into account in this paper. 
Because there are much more terms in the 4-loop order 
in Fig. 4 and they are almost equivalent (see §4.4), the 
unimportance of the non-Gaussian term is remarkable in 
the 4-loop order. This is partly owing to the fact that the 
4-th order vertex T — £fc|G(fc)| 4 is significantly reduced 
by the self-energy arising from the spin fluctuation as 
well as the SC fluctuation. Although the role of higher 
order non-Gaussian terms is still unclear, we believe that 
they do not seriously affect the qualitative conclusions 
obtained in this paper. 
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